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Abstract. We constructed hydrodynamical model atmospheres for mid M-type main-, as well as pre-main-sequence 
(PMS) objects. Despite the complex chemistry encountered in these cool atmospheres a reasonably accurate 
representation of the radiative transfer is possible, even in the context of time-dependent and three-dimensional 
models. The models provide detailed information about the morphology of M-type granulation and statistical 
properties of the convective surface flows. In particular, we determined the efficiency of the convective energy 
transport, and the efhciency of mixing by convective overshoot. The convective transport efficiency was expressed 
in terms of a n equ ivalent mixing-length parameter q:mlt in the formulation of mixing-length theory (MLT) given 
bv iMihala j il978ll . CfMLT amounts to values around ~ 2 for matching the entropy of the deep, adiabatically 
stratified regions of the convective envelope, and lies between 2.5 and 3.0 for matching the thermal structure 
of the deep photosphere. For current spectral analysis of PMS objects this implies that MLT models based on 
Qmlt = 2.0 overestimate the effective temperature by 100 K and surface gravities by 0.25 dex. The average thermal 
structure of the formally convectively stable layers is little affected by convective overshoot and wave heating, 
i.e., stays close to radiative equilibrium conditions. Our models suggest that the rate of mixing by convective 
overshoot declines exponentially with geometrical distance to the Schwarzschild stability boundary. It increases 
at given effective temperature with decreasing gravitational acceleration. 
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1. Introduction 

The increasing number of stars, brown dwarfs, and ex- 
trasolar planets of spectral class M or later discovered by 
infrared surveys and radial velocity searches has spawned 
a great deal of interest in the atmospheric physics of these 
objects. Their atmospheres are substantially cooler than 
the solar atmosphere, allowing the formation of molecules, 
or even liquid and solid condensates. Convection is a 
ubiquitous phenomenon in these atmospheres shaping 
their thermal structure and the distribution of chemical 
species. Hydrodynamical simulations of solar and stel- 
lar granulation including a realistic description of ra- 
diative transfer have become an increasingly powerful 
and handy instrument for studying the influence of con- 
vective flows on the structure of late-type stellar at- 
tend offprint requests to: Hans-Giinter Ludwig 



mosp heres as well as on the formation o f their spectra 
(e.g. iNordlundl Il982t ISteffen et alJll989t Fchan fc Sofial 
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Vogler fc SchiissleJl200:j:lR"obinson et al.ll2004l) . Here we 



report on efforts to construct hydrodynamical model at- 
mospheres for mid M-type objects. The spectral type just 
borders the temperature where the formation of conden- 
sates becomes important. The motivation of this investi- 
gation was twofold: first, pre-main-sequence (PMS) evolu- 
tionary models of M-type star s and br own dwarfs based on 
mixing -length theory (MLT. lBoh~Vitense 195^) to de- 
scribe the convective energy transport depend sensitively 
on the poorly constrained mixing-length parameter^ Q?mlt 

^ The ratio between the mixing-length and local pressure 
scale height. 
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l)Baraffe et al.ll2002|) . Our hydrodynamical models repre- 
sent convection essentially from first principles, and are 
free of the uncertainties of MLT allowing to put the 
stellar models on a firmer footing. Second, the distribu- 
tion of dust clouds in cool brown dwarfs depends on the 
efficiency of mixing of their atmospheres by convective 
overshoot llAckerman fc MarlevI lioOlt lAllard et aLll2003t 
iHeUing et all 120041) . For other work on the modeling of 
dust cloud for mation in very low mass stars and b rown 
dwarfs, see e.g. lCooper et al.l l|200,'^ and lTsuiil l|20n,'J) . and 
references therein. While by construction MLT cannot de- 
scribe convective overshoot it is naturally represented in 
our three-dimensional hydrodynamical models. 

The present investigation is extending a previous study 
of an M-dwarf atmosphere by iLiidwig et all l|2002l here- 
after LAH) to PMS objects at lower surfac e gravity . A pre- 
liminary account of the results was given in iLudwid 12003). 
We start in Sect.|2]with an overview of the model construc- 
tion, in particular related to approximations we adopted 
in the radiative transfer. In Sect.|3|we discuss the general 
morphology of the convective flows in the M-type atmo- 
spheres and present some statistical properties. In Sect.0] 
we provide estimates of the efficiency of the convective 
energy transport in terms of an effective mixing-length 
parameter, and discuss consequences for the analysis of 
PMS M-type objects in the framework of present standard 
model atmospheres. We continue in Sect.[Slby characteriz- 
ing the properties of the atmospheric mixing found in the 
hydrodynamical models, and conclude with final remarks 
in Sect. El In our investigation we take repeatedly recourse 
to the solar atmosphere as standard benchmark. 



2. Model overview 

Figure H illustrates the positions of our hydrodynamical 
model atmospheres in the Teff-logg-plane: three M-type 
models are located close to T^s = 2800 K with log g=3.0, 
4.0, and 5.0 which form a logg-sequence. In the follow- 
ing we shall refer to them as models C3, C4, and C5, 
respectively. To assess temperature effects, a 500 K hot- 
ter M-type model was constructed at Toff =3280 K and 
log (7=4.0. For further comparison we also considered a 
solar (in terms of its hydrodynamical properties) model S 
at Toff =5640 K, log.g=4.44, and a model SG of a subgiant 
at 4610 K and log 5=2.94. 

For investigating the influence of the position of the 
upper boundary condition we constructed an additional 
model HX with the same atmospheric parameters as 
model H4 but extending 270 km (corresponding to 3.7 Hp) 
higher up than model H4. The flow in the extended region 
exhibits larger fluctuations than encountered in deeper 
layers favoring the formation of sharp flow features. For 
reasons of numerical stability we had to increase the nu- 
merical viscosity so that the model is not fully differen- 
tially comparable to model 114. Nevertheless, it should 
give an indication of the level of the influence of the up- 
per boundary, in particular when the upper boundary is 
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Fig. 1. The radiation- hydrodynamics models in the effec- 
tive temperature-gravity plane (cubes). The models are 
labeled by their IDs used for reference in this paper. The 
approximate spectral class is indicated at the temperature 
axis. 



located close to the convectively unstable region. In plots^ 
that follow, we depict model HX always as triple-dot- 
dashed line without labeling it by its name like the other 
models. Naturally, the behavior of model HX closely fol- 
lows model H4 in the deeper layers so that its connection 
to model H4 is readily apparent. 

All models were evolved until a thermally and dynam- 
ically relaxed state was reached. All models except HX 
have 125x125x82 grid points (XxYxZ direction), HX has 
125x125x102 grid points due its larger vertical extent. The 
numerical grid is equidistant in x- and y- direction while 
in vertical z-direction the grid spacing is to first order cho- 
sen to provide the same number of grid points per pressure 
scale height. In addition, the resolution is increased in lay- 
ers around continuum optical depth unity if a steep verti- 
cal temperature gradient is present. All models have solar 
chemical composition. Table ^ summarizes their proper- 
ties. 

The overall methodology applied in this work is the 
same as in LAH, and we refer the reader to this paper for 
details beyond the short description we provide here. 

The radiation-hydrodynamics (RHD) simulations 
were performed with a convect ion code developed by 
A. Nordlund and. R.F. Stein fsee lStein fc Nordlundll99gL 
and references therein). The code solves the hydrodynam- 
ical equations of compressible gas dynamics coupled with 
non-local radiative transfer in three spatial dimensions. 
The time-independent radiative transfer is treated assum- 
ing strict LTE. The wavelength dependence of the radia- 
tion field is represented by a small number of wavelength 
bins. Open lower and upper boundaries, as well as periodic 
lateral boundaries are assumed. The effective temperature 
of a model (i.e., the average emergent radiative flux) is 
controlled indirectly by prescribing the entropy of inflow- 
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Table 1. The RHD models discussed in the paper: ID is the identifier used to refer to a model in this paper, Teff 
the effective temperature of the model including an estimate of its RMS fluctuations, logg the preset gravitational 
acceleration, Size is the geometrical size of the computational domain (XxYxZ, where Z denotes the vertical, X and 
Y the horizontal directions), H^'^^ the pressure scale height at the surface, Pgurf the pressure at the surface, the 
maximum vertical RMS velocity in the convective layers, (5/rms/^ the relative intensity contrast, omlt (evo) the mixing- 
length necessary to match the asymptotic entropy (see Sect.^, q;mlt (phot) the mixing-length parameter necessary 
to match the model's temperature in the deep photosphere, iJfex the scale height of the decline of the atmospheric 
mixing rate (see Sect. 01, and Modelcode an internal model identifier. Parenthesis indicate uncertain values. 



ID Toff logg Size i/f,"'^' P^urf vf^^ 5/rms// Qmlt Qmlt -fffox Modelcode 

[K] [Mm] [Mm] logio [ms-^] [%] (evo) (phot) [Hp] 
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ing material at the lower boundary. Magnetic fields and ro- 
tation are neglected. Opacities and equation-of-state have 
been adapted to the conditions encountered in M-type at- 
mospheres. The equation of state includes the ionization of 
H and He , as well as H2 molecule formation according to 
Saha-Boltzmann statistics. II2 molecule formation is the 
thermodynamically most important process in the M-type 
atmospheres. The opacities include contributions of molec- 
ular lines but neglect contributions of dust grains which 
is a good approximation at the temperatures prevailing in 
the models. The opacities were extracted from the opacity 
data base of the PHOENIX model atmosphere code (for a 
description of PHOENIX and corresponding opacities see 
iHauschildt et al.lll999tlFerguson et al.ll2005|) . 

2.1. Radiative transfer 

We want to derive quantitative estimates of the mixing 
by convective overshoot, as well as obtain a measure of 
the efficiency of the convective energy transport. For ad- 
dressing these issues, the RHD models have to give a 
reasonably accurate representation of the actual atmo- 
spheric conditions. Here we are particularly concerned 
about the radiative energy transport, which is compli- 
cated by the huge number of molecular absorption lines. In 
our RHD models, we use a multigroup technique (dubbed 
Opacity Binning Method, hereafter OEM) for modeling 
the radiative energy exchange which employs four groups 
for represe nting the wavelength dependeii c e of the radi- 
ation field llNord lund"l982': Lvidwid Il993 iLudwig et aP 
119941 IVogler et al. 2004. LAH). The wavelength groups 
have been optimized for an atmosphere at rcff=2900K 
and log 5=5.3. Figure [3 illustrates the accuracy which is 
achieved with the OEM for the present models. We com- 
pare ID MLT model structures (aMLT=l-0) in radiative- 
convective equilibrium computed with the OEM approxi- 
mation and high-precision opacity sampling. While there 



are differences between the atmospheric structures, the 
OEM nevertheless provides a significant improvement 
with respect to a simple grey approximation. Temperature 
differences get larger as one moves away from the at- 
mospheric parameters the OEM was optimized for, and 
reach up to 250 K in the model at 2800 K and log 5=3.0. 
However, in the present context it is not so much the ab- 
solute temperature error as the change of the temperature 
gradient which is relevant. 

Convection is driven by buoyancy forces whose dynam- 
ical effects scale with the entropy gradient. The thermal 
structure is controlled by a balance between radiative and 
advective (due to compression or expansion of mass ele- 
ments) heating (or cooling) of which the Peclet number 
(see appendix ^ for its definition and computation) is 
a convenient dimensionless measure. In the deeper pho- 
tospheric layers, the OEM profiles shown in Fig. |2 have 
steeper temperature gradients than the profiles based on 
opacity sampling. For a given velocity, this makes the rate 
of temperature change of a mass element moving in ver- 
tical direction larger. In other words, the time scale of 
advection related temperature changes becomes shorter, 
and the Peclet number larger as long as the radiative time 
scales remain the same. The opposite behavior is present 
in the higher photospheric layers. The typical Peclet num- 
bers turn out to be by a factor of 1.2 times larger for the 
log g=5.0, and by up to a factor of 2 times smaller for the 
log (7=3.0 model comparing OEM to the opacity sampling 
stratifications. 

To mitigate the effects of the shortcomings, primarily 
related to the OEM approximation in the radiative trans- 
fer, we took a differential approach when measuring model 
properties: whenever possible we compared RHD and hy- 
drostatic model atmospheres based on the same OEM ra- 
diative transfer scheme since we were interested to study 
the systematic change of model properties with effec- 
tive temperature and gravitational acceleration. However, 
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Fig. 2. Comparison of the pressure-temperature structure 
between ID hydrostatic model atmospheres in radiative- 
convective equilibrium based on the OBM (solid lines) 
with four wavelength groups, and PHOENIX models 
based on direct opacity sampling employing several 10^ 
wavelength points (dash-dotted) . Shown are two examples 
at Teff=2800K, log(jf=3.0 and 5.0, respectively. For clar- 
ity, the log (7=3.0 models were shifted by -1-1000 K. Also 
shown is a model employing grey radiative transfer (dot- 
ted line). The pressure is given in units of the pressure at 
Rosseland optical depth unity Psurf- 

there remains the possibility that even this differential ap- 
proach cannot fully eliminate the impact of the artificial 
shift of the Peclet number introduced by the OBM ap- 
proximation among the models. This has to be kept in 
mind when interpreting results later. 

3. Morphology of granulation in M-type objects 

FigureElshows an inter-comparison of the granulation pat- 
terns typically encountered during the temporal evolution 
of our RHD models. The first thing to note is that sur- 
face convection in M-type objects produces a granular 
pattern qualitatively resembling solar-type granulation: 
bright extended regions of up-welling material which are 
surrounded by dark concentrated lanes of down-flowing 
material. The dark lanes form an interconnected net- 
work. Looking more closely, however, granules are less reg- 
ularly delineated in M-type objects. The inter-granular 
lanes show a higher degree of variability in terms of 
their strength - in particular in comparison to the so- 
lar model S, to lesser extend in comparison to the sub- 
giant model SG. A feature which is uncommon in the solar 
granulation pattern are the dark "knots" found in or at- 
tached to the inter-granular lanes prominent in the M-type 
models C5 and C4. The knots are associated with strong 
downdrafts which carry a significant vertical component 
of angular momentum. 

We have no convincing explanation at hand why con- 
vective flows in M-type atmospheres tend to form such 
vortical structures. It is unlikely that they are (as sus- 
pected by the referee) an artifact related to the boundary 



condition since the structures do not reach up into lay- 
ers close to it. Moreover, the test model HX shows also 
the knots despite the fact that its upper boundary is lo- 
cated about 10 Hp above the continuum forming layers. 
The presence or absence of knots may rather be related 
to the level of horizontal shearing in the layers around 
optical depth unity where convective driving is usually 
strongest and is important for the formation of flow struc- 
tures. Figured illustrates that under solar-like conditions 
(models S and SG) shear flows are common in the sur- 
face layers, and are much more pronounced than in the 
M-type atmospheres. Such shear flows make the forma- 
tion of structures extending in vertical direction difficult, 
and may be the reason why the knots do not appear in 
the solar-like models and are less developed in models H4 
and C3. 

As a side-point we would like to remark that it is 
not immediately clear why the knots appear in fact dark. 
However, centrifugal forces in the vortices are evidently 
not strong enough to lead to an evacuation of their interior 
on a level that would allow radiation from deeper, hotter 
layers to escape and let them appear brighter than their 
surroundings. The conditions are not as extreme as for 
bright flux tubes in the solar photosphere where magnetic 
pressure allows for a substantial degree of evacuation. 

The width of the inter-granular lanes relative to the 
typical granular size is smaller in our M-type objects of 
higher gravity. Inspecting the velocity field (not shown) 
in vicinity of the continuum forming layers shows less 
pronounced size differences. This indicates that the rel- 
atively broader lanes in the solar case are the result of 
a stronger smoothing of the temperature field due to a 
more intense radiative energy exchange, i.e., the effectively 
smaller Peclet number of the flow around optical depth 
unity in the hotter objects. Figure ^ shows an overview 
of the mean T-r (the temperature was averaged over time 
and surfaces of equal optical depth) relations found in the 
RHD models. One recognizes that in the M-type objects 
the thermal structure is influenced by convection to much 
lower optical depths than in the solar- type stars. 

The temperature model HX is on the scale of the plot 
identical to model H4 in deeper layers but shows a notice- 
able deviation at low optical depth close to the upper limit 
over which the stratifications have been averaged. While 
perhaps not surprising considering the different placement 
of its upper boundary, the difference may be in part traced 
back to the different viscosity which damps the velocity 
field and leads to less convective heating in the atmo- 
spheric layers which are convectively unstable (see Sect.0] 
for a discussion of the interplay of convection and radia- 
tion in the atmospheric layers). 

3.1. Horizontal scales 

Primarily due to the variation of the gravitational acceler- 
ation (by a factor of 100) the convective cells in our models 
span a substantial range in geometrical size. However, tak- 
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Fig. 3. Granulation patterns of the six RHD models. The model IDs are given in the upper left corners (see Tab.^l of 
the images. Shown are snapshots of the emergent intensity in the continuum. For each image a separate grey scale is 
used, lighter shades correspond to higher intensities. The relative intensity contrast of the particular image is stated, 
and the bar indicates a length of 10 x H™'^^ . 



ing the pressure scale height at optical depth unity _ffp"''^ 
as reference the variation is largely reduced. The bars in 
Fig. 13 are placed to suggest that the horizontal size of the 
cells indeed roughly scales with iJp . This is quantified in 
Fig.Elfor the intensity pattern and Fig.|B|for the pattern of 
the vertical velocity. The maximum power in the spectra 
of the intensity pattern of the M-type objects lies between 
5 and 8 iJp"''^, the maximum power of the velocity pattern 
between 3 and 6i?p"''^ In both cases the spectra of the 
solar- type objects are slightly but noticeable shifted to- 
wards smaller wavenumbers (larger spatial scales). While 
the relation between M-type and solar-type objects is the 
same in both diagnostic variables, it came as a surprise - 
at least to the authors - that both variables do not pro- 
vide the same value for the typical cell size. The different 
slopes in the intensity and velocity spectra towards larger 
scales might be related to this finding. If one considers a 
pure random pattern of a given characteristic scale, a slope 
of unity (in the chosen representation of power) is to be 
expected towards larger scales. In this case the signal at 
large scales is the result of a mere random superposition 
of residual contributions from smaller scale features. The 
intensity spectra follow this random model quite closely 



while the velocity spectra show noticeably larger devia- 
tions with steeper slopes. We do not have an explanation 
at hand. However, considering the range of stellar param- 
eters covered by our models, the shape of spectra show a 
large degree of similarity. In particular, the typical gran- 
ular scales (in intensity) turn out to be the same within a 
factor of two lying between 5 and 10 ifp"''*for all objects. 

3.2. Velocities and turbulent pressure 

Figures H and |H1 show the run of the root-mean-square 
(RMS) vertical as well as horizontal velocity component, 
and of the turbulent pressure, respectively. The averages 
were taken over time and fixed geometrical height. As is 
evident from the figures, the vertical velocity and tur- 
bulent pressure follow each other rather closely. In the 
higher atmospheric layers, the velocity field - especially 
in the models with more vigorous convection - is domi- 
nated by horizontal motions. Test model HX closely fol- 
lows model H4 in the deeper layers but shows reduced am- 
plitudes in the higher layers. As indicated previously we 
interpret this as primarily a consequence of the increased 
numerical viscosity in model HX, and not so much due to 
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Fig. 4. Mean temperature (averaged over time and sur- 
faces of constant optical depth) relative to effective tem- 
perature for the RHD models as function of Rosseland 
optical depth. For clarity the curves are shifted by 
{1, 2, 3, 4, 5} X 0.05 for the models C4, C3, H4/HX, S, and 
SG, respectively. 
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Fig. 6. Spatial power spectra of the vertical velocity as a 
function of horizontal wavenumber k. The velocities are 
taken from the layer where the convcctive velocities reach 
their maxima. The curves are labeled with the IDs of the 
models (see Tab. the curve labeled "1" indicates a 
power-law with a slope of one. The grey stripe indicates 
the range of horizontal scales within which the power max- 
ima of the M-type models are located. 
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Fig. 7. Root-mean-square velocities for the RHD models 
as a function of gas pressure relative to the individual 
surface pressures Psuri (see Tab.^. Black lines depict the 
vertical, grey lines the horizontal velocity component. 



the larger extend of the model. Generally, the velocities 
are smaller in the M-type objects then in the solar-like 
objects. Due to their substantially lower Toff in compari- 
son to the solar-type models, the requirement to transport 



the nominal energy flux is already met by convection at 
lower velocities in the M-type models, leading to the over- 
all lower "hydrodynamic activity" in these objects. 
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Fig. 8. Turbulent pressure {pv^Vz) as a function of gas 
pressure relative to the individual surface pressures. The 
turbulent pressure quite closely follows the vertical veloc- 
ity as depicted in Fig. |7| 



3.3. Horizontal fluctuations 

Figures El and EH show the relative spatial and temporal 
RMS fluctuations of temperature and pressure at given 
geometrical height. In the M-type models the tempera- 
ture fluctuations stay at a very modest level nowhere ex- 
ceeding 4 % - even including model HX. In the optically 
thin layers, model C4 shows systematically larger temper- 
ature fluctuations than model H4 of higher T^s but of the 
same surface gravity. This is an imprint of the reduced 
capacity of the radiation field in model C4 to smooth hor- 
izontal temperature differences. The pressure fluctuations 
reach larger values than the temperature fluctuations, and 
show a systematic increase with height. Note that in our 
dwarf model C5, the pressure fluctuations only reach a 
very modest level of about 4%. From this we expect that 
in cooler, dust-forming main-sequence objects, thermody- 
namic fluctuations are even smaller so that dust formation 
conditions vary little in a given layer. Model HX shows a 
rapid increase of fluctuations with height which is typically 
found when the flow field is dominated by wave motions. 
The correspondence to model H4 in the overlapping region 
is quite good indicating that fluctuations in temperature 
and pressure are little affected by the location of the upper 
boundary. 

3.4. Spatial correlation of vertical velocity and entropy 

In this section, we want to give an overview of some two- 
point correlations found in our models. The width of two- 
point correlations in vertical direction has been considered 
as measure of the mean-free-path of mass elements enter- 
ing MLT ~ the mixing-length A, and as such has been 
the target of many investigations. The (linear) correlation 
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Fig. 9. Relative RMS temperature fluctuations in hori- 
zontal planes as a function of gas pressure relative to the 
individual surface pressures Psmf (see Tab.P). The clipped 
maxima in the curves for models S and SG reach 0.23 and 
0.30, respectively. 
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Fig. 10. Relative RMS fluctuations in horizontal planes as 
a function of gas pressure relative to the individual surface 
pressures. 

coefficient of a quantity x between two layers located at 
heights zi and Z2 is given by 



C [xi,X2] 



{X1X2) - {Xi) (X2) 



(1) 



where a^. is the standard deviation of x at height z^, and 
the angular brackets denote the average over time and hor- 
izontal position. In a seminal paper, Chan fc Sofia (19i3) 
found in hydrodynamical models of stratified efficient con- 
vection that the fuU-width-at-half-maximum (FWHM) of 
the correlation function of the vertical velocity, as well as 
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temperature, scaled with the pressure scale height and not 
densit y scale height. In a follow-up study, IChan fc Sofial 
l)l989(l showed that the result is robust ag ainst variations 
of the ratio of the specific heats 7 of the gas. lSingh fc ChanI 
1)1993(1 found that the width o f correlation chan ges mod- 
erately with Prandtl number. iKim et al ] lll995^ studied 
the case of inefficient convection in the Sun. Their model 
included radiative transfer effects (in diffusion approxima- 
tion) as well as effects of the ionization of hydrogen and 
helium. Under the conditions studied by Kim and collab- 
orators, the width of the correlation function of the veloc- 
ity scaled with pressure and density scale height, while the 
width of th e temperature correlation with neither of both. 
iKim et all ha d to restrict the i r investigati on to optically 
thick regions. iRobinson et"al] l|2003l l2004(l included also 
optically thin layers in their hydrodynamical models of 
the Sun and a subgiant star (at rcff=4990K, log (7=3.37). 
The radiative transfer was treated in grey Eddington ap- 
proximation. Robinson and collaborators found a complex 
height-dependence of the width of the correlation of the 
vertical velocity and entropy fluctuations. Especially the 
last result shows that we cannot expect to find a universal 
behavior of the correlations in stratified convection under 
all possible circumstances - in particular, in the case of 
inefficient convection which is most important for stellar 
structure models. However, there might be still hope to 
find general trends which may serve as buildings blocks 
for an improved treatment of convection beyond MLT. 

Figures [TTl and [T^ provide an overview of the velocity 
and entropy correlations found in our models. Our mod- 
els contribute to the ongoing discussion in a twofold way: 
they cover a large range of stellar parameters and include 
an elaborate treatment of the radiative transfer in the op- 
tically thin layers. The overall behavior of the correlations 
is complex, but some general features can be identified: i) 
considering that the models cover about a factor of two 
in Teff and one hundred in logg, a certain uniformity in 
the width of the distributions in the convection dominated 
layers is apparent; ii) in the sub-photospheric layers the 
entropy correlation is more peaked than the velocity cor- 
relation; iii) with the exception of model C5 the velocity 
correlation become broader in radiation dominated layers; 
iv) the widths of the correlations tend to shrink with in- 
creasing depth, possibly towards an asymptotic limit (dif- 
fe rent for velocity and entropy) - in line with the findings 
of lChan S(^ lll 98flh : v) the width of the entropy corre- 
lation of the solar model S and subgiant SG passes through 
a pronounced minimum around optical depth unity; this 
is accompanied by anti-correlations signifying the thermal 
behavior of overshooting motions; vi) the degree of unifor- 
mity of the correlations is not improved if plotted on the 
density scale. 

The limiting width of the correlation functions is fairly 
well defined and given in the panels of Figs.lTTIandlT^ The 
width of the velocity correlation lies in the same range as 
values of the mixing-length parameter associated with cer- 
tain model features which we are going to discuss later. 
However, one must keep in mind that the width of the 



correlations varies substantially - sometimes even dramat- 
ically - in the models, and that the widths are different for 
different quantities. In our opinion, a direct association of 
a width of a correlation with a mixing-length parameter 
is an over-simplification of the actual situation. 

4. Convective energy transport 

Convection is an important energy transport mechanism 
in M-type stars. In standard model atmospheres it is 
treated in the framework of MLT. In this section, we 
want to address the question whether the simplistic MLT 
is actually capable to provide a sufficiently accurate de- 
scription of the convective energy transport under con- 
ditions encountered in M-type atmospheres. Figure 1131 
shows a comparison of the entropy structure of the RHD 
model atmospheres and standard ID hydrostatic mod- 
els in radiative-convective equilibrium assuming different 
mixing-length parameters. Figure IT^ depicts correspond- 
ing temperature profiles which allow to approximately 
translate entropy differences among the profiles in Fig. 1131 
to temperature differences. 

To calculate the entropy profiles of the RHD mod- 
els, they have been averaged temporally and horizontally 
on surfaces of constant optical depth. This procedure en- 
sures a particularly good preservatio n of the energy trans - 
port properties of the RHD models l|Steffen et al.lll995|) . 
Moreover, it reduces the "smearing" of vertical gradients 
by plane-parallel oscillations which occurs when averaging 
over fixed horizontal planes. Besides the mixing-length pa- 
rameter itself, MLT contains a number of further "hidden" 
parameters intrinsic to the specific formulation of MLT 
which was chosen. We emphasize that a well-defined cal- 
ibration of the mixing-length parameter must always be 
given with reference to the specific formulation i n opera- 
tion. H ere, we are using th e formulation given bv lMihalasI 
l|l978l) . See lLudwig et alJ fl999 ) for details of the imple- 
mentation. 

As already remarked earlier, Fig.ll3lillustrates that the 
sensitivity of the structure of the standard models to the 
mixing-length parameter increases with decreasing grav- 
ity as well as increasing Teff. Model C3 shows a sensitivity 
of the entropy in the deep, adiabatic layers which is com- 
parable to the sensitivity of solar MLT models. Figure [T^ 
also illustrates that the convectively unstable layers (with 
entropy gradient ^ > 0) extend to small - spectroscopi- 
cally important - optical depths. M-type atmospheres of- 
fer the opportunity to study convection under optically 
thin conditions. 

The general role of convection can be described as fol- 
lows; in the convectively unstable layers - here comprising 
also parts of the optically thin layers - the thermal struc- 
ture is the result of a competition between adiabatic heating 
and radiative cooling because a temperature structure in 
adiabatic equilibrium would be hotter than in radiative 
equilibrium. The mixing-length model (with omlt — 2.0) 
depicted by the dashed-dotted line for case H4 is intended 
to illustrate this (see Figs. El ^2 and[T7|l: in this model 
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Fig. 11. Contour plots of the two-point spatial correlation 
functions of the vertical velocity (left panels), and entropy 
(right panels) for the models C3, C4, and C5. The pres- 
sure relative to the surface pressure of the two involved 
height levels (cf. Eq. (Q)) are given on abscissa and or- 
dinate. The height coordinates are interchangeable due 
to the symmetry of the correlation function. The number 
in the lower right corner of each panel gives the FWHM 
of the correlation (in units of the local pressure scale 
height i/p) in vicinity of the lower boundary of the com- 
putational domain. Contour lines are given for values of 
{-0.5, -0.25, -0.125, 0.125, 0.25, 0.5}. Lines with positive 
contour values are depicted by solid lines, negative ones 
by dashed lines. Dotted lines are for orientation and have 
a distance of Aln(P/Psurf) — 1- The central dotted line 
coincides with the maximum of the correlation function 
which is normalized to one. 



we artificially switched off the convective motions in the 
layers with logTross < ^1- This suppresses the convective 
heating and forces the temperature to adjust to radiative 
equilibrium conditions. As evident from the figures, this 
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Fig. 12. Same as Fig. HH for models SG, S, and H4. 



leads to a substantial drop of entropy and temperature. 
In the convectively stable layers (with entropy gradient 
^ < 0) the situation is reversed, and the temperature is 
controlled by a balance between adiabatic cooling and ra- 
diative heating - as far as the RHD models are concerned. 
The situation is different for the MLT models where by 
construction no convective (overshooting) motions take 
place in the formally stable layers, and the temperature 
is determined by the condition of radiative equilibrium 
alone. 

For decreasing log 5 at given Toff, the models tend to 
stay closer to radiative equilibrium conditions in the opti- 
cally thin layers. Two factors reduce the efficiency of the 
convective energy transport: first. Fig. 1151 shows that the 
opacity does not vary radically among the M-type mod- 
els. Hence, lower gravity models exhibit lower densities at 
given optical depth, which reduces the thermal energy that 
can be transported per unit volume, rendering the convec- 
tive transport of heat more difficult. Second, the pressure 
scale height increases with decreasing gravity while typi- 
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Fig. 13. Entropy as a function of Rosseland optical depth 
of the RHD models (thick solid lines) in comparison 
to standard mixing-length models (thin solid lines). For 
each RHD model three MLT models are plotted with 
aMLT=l-5, 2.0, and 2.5. The entropy of the MLT models 
behaves monotonically with auLT, the q;mlt=1-5 model 
having the lowest entropy in the optically thin, and the 
highest in the optically thick layers. Model H4 has been 
offset by +0.2 entropy units for clarity, and one special 
MLT model with q;mlt=2.0 has been added (dashed- 
dotted line; see text). The dashed lines depict the value 
of entropy present in the adiabatically stratified regions of 
the convective envelope. Numbers indicate mixing-length 
parameters necessary to match the RHD structure by 
MLT models. 



cal convective velocities increase only modestly. This in- 
creases the time scale over which vertically traveling mass 
elements change their temperature due to adiabatic ex- 
pansion or compression. The radiative time scale in the 
optically thin regions, on the other hand, is independent 
of the spatial scales and mass density which again leads 
to a shift of the thermal balance towards radiative equi- 
librium conditions. 

As evident from Fie. 1161 the pressure-temperature de- 
pendence of the adiabatic gradient counter-acts the trend 
towards radiative equilibrium conditions in our M-type 
models of lower gravity. In models C3 and H4, the adia- 
batic gradient is close to its minimum favoring convection 
due to the formation of H2 molecules. Nevertheless, the 
appreciable sensitivity of the MLT models to the mixing- 
length parameter at lower gravities indicates that convec- 
tion and radiation operate with comparable efficiency. 

In Figs. El and El we compare groups of models of 
the same Toff and logg, and consequently the models 
in each group are constrained to similar temperatures 
and entropies in vicinity of optical depth unity. Since 
the entropy- as well as temperature-gradient at this lo- 
cation depend on the mixing-length parameter, this leads 
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Fig. 14. Temperature as a function of Rosseland optical 
depth of the RHD models (thick solid lines) in compar- 
ison to standard mixing-length models (thin solid lines). 
The MLT models are calculated assuming aMLT=l-5. For 
model H4 a special MLT model with aMLT=2.0 has been 
added (dashed-dotted line; see text). For clarity, models 
C4, C3, and H4 have been offset by 500, 1000, and 1500 K, 
respectively. 
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Fig. 15. Average pressure-temperature profiles of the 
RHD models (thick lines) overlayed on contours (thin 
lines) of the (logj^p) Rosseland mean opacity. 



to the crossing of the various profiles around Tj-oss = 1- 
Deviations from this behavior - in particular shown by the 
RHD models - come about by residual differential changes 
of the opacities and thermodynamic properties of the mat- 
ter among the models. In the case of the RHD models, 
horizontal fluctuations together with non-linearities in the 
material functions add to the deviations. 

Figure El shows that irrespective of the choice of the 
mixing- length parameter ckmlt, no MLT model is capable 
to match the whole average thermal profile of a hydrody- 
namical model - at least within the framework of the MLT 
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Fig. 16. Average pressure-temperature profiles of the 
RHD models (thick lines) overlayed on contours (thin 
lines) of the adiabatic gradient. 



formulation adopted here. We nevertheless give estimates 
of the mixing-length parameter matching certain features 
of the hydrodynamical structure. Multi-dimensional hy- 
drodynamical convection models provide the value of the 
entropy asymptotically reached in the deep, adiabatically 
stratified layers of the convective envelope, since it is iden- 
tical to the entropy of upflowing material in the deepest 
regions of the hydrodynamical model (for a d i scuss ion of 
the scenario underly ing: this notion see lStefFenl l(l99.ll) . and 
iT.udwig et all l)l999t) 1. We find a mixing-length parame- 
ter of 1.5, 2.1, 2.1, and 1.85 for models C5, C4, C3, and 
H4, respectively, to match the asymptotic entropy. This 
Q^MLT is most relevant for stellar structure models since 
the asymptotic entropy has a large influence on the radius 
of a convective stellar object. The value for model C5 (1.5) 
is highly uncertain since the sensitivity to qmlt is very 
small ~ however, its precise value is also not particularly 
important since the asymptotic adiabat hardly depends 
on aMLT- Attributing little weight to model C5, we find 
that the typical q:mlt suitable for evolutionary models 
lies around « 2 for mid to late M-type atmospheres. This 
value is not very different from the value found for the 
Sun (« 1.8), but we do not consider this as an indication 
of an universal value for q;mlt a-s was already discussed in 
LAH. 

Now we turn to the optically thin parts of the atmo- 
sphere relevant for spectroscopy. We find that mixing- 
length parameters for matching the temperature in the 
range —2 < logTross < — 1 of 2.5, 2.5, 2.8, and 3.0, for mod- 
els C5, C4, C3, and H4, respectively. This range in optical 
depth was chosen as representative of the deeper photo- 
sphere. However, the match in this part does not imply 
a match over the whole optically thin region. Again, the 
value for model C5 is uncertain, but not particularly im- 
portant. All in all we obtain a range of omlt — 2.5 . . . 3.0 
when matching the thermal profile of the RHD models. 
The temperature structure of the hydrodynamical models 
in the uppermost, convectively stable part of the atmo- 



sphere does not deviate much from radiative equilibrium 
profiles judging by extrapolating from the available MLT 
models. This indicates that cooling by convective over- 
shoot or heating by waves is not very efficient in the layers 
immediately adjacent to the surface convective zone. 

Model H4 is an exception from the general trend since 
it becomes convectively stable much earlier than one might 
expect from extrapolating the MLT models. The excep- 
tional behavior is due to the fact that in this model the 
adiabatic gradient changes very little along the profile in 
the upper photosphere (see Fig. The same holds for 
the actual temperature gradient in the convectively un- 
stable part, making the exact location of the transition 
from convective instability to stability very sensitive to 
the actual run of the temperature. Here, we have an ex- 
ample where second order effects can enhance differences 
to MLT models. Nonetheless, overshooting and wave heat- 
ing, again, show little impact on the temperature gradient 
in the (not very extended) stable zone. 

A comparison of models H4 and HX in Fig. El re- 
veals that the modest temperature differences in the op- 
tically thin layers displayed in Fig. E| correspond to a 
sizable entropy difference which corresponds to a decrease 
of the mixing- length parameter of about 0.5 relative to 
model H4. Such a change would bring the atmospheric 
value of the mixing-length parameter of model H4 closer 
to the other M-type models. As argued before the change 
in the thermal structure of model HX can only be par- 
tially attributed to a systematic infiuence of the upper 
boundary condition. Nevertheless, taking the change as 
an estimate of the uncertainty of the derived atmospheric 
mixing-length parameter, one would still find a rather 
high value of the mixing-length parameter for model H4 
of at least 2.5. The mixing-parameter for matching the 
asymptotic entropy remains the same between models H4 
and HX. 

In LAH we discussed the mixing-length parameter nec- 
essary to match the maximum vertical RMS velocity of the 
hydrodynamical model C5 and found a value of 3.5. While 
we do not perform a similar matching here, Fig. llTl clearlv 
shows that for all RHD models a value substantially larger 
than 2.5 - including model HX - is necessary to match the 
maximum velocity. 

Comparing the overall situation that we encounter in 
M-type objects to stars of roughly solar effective temper- 
ature (and of solar composition) wc find that the tran- 
sition from convectively to radiatively dominated energy 
transport happens more gradually in M-type objects. This 
is ultimately linked to the different temperature sensitiv- 
ity of the dominant opacity (H and H~ bound-free and 
free- free absorption versus TiO and H2O molecular line 
plus H^ bound- free and free-free absorption), and domi- 
nant thermodynamic process (H recombination versus H2 
molecular formation) encountered in the two regimes of ef- 
fective temperature. The values of omlt we obtained here 
when matching the asymptotic entropy are not too differ- 
ent from the values obtained for solar- type stars. However, 
in solar-type stars a calibrated MLT model merely pro- 



12 



Ludwig, AUard, Hauschildt: Radiation-hydrodynamics of M-type atmospheres 



30001 — r 



2500 




log 10 '^To 



Fig. 17. RMS vertical velocity component of the RHD 
models (thick solid lines) in comparison to convective ve- 
locities from MLT models with aMLT=l-5, 2.0, and 2.5 
(thin solid lines) as a function of optical depth. The veloc- 
ities from MLT increase in the optically thin layers mono- 
tonically with increasing cimlt- For model H4 a special 
MLT model with aMLT=2.0 has been added (dash-dotted 
line; see text). Models C4, C3, and H4 have been shifted 
by 500, 1000, and 2000 ms~^, respectively. The downturn 
of the velocity in models C4, C3, and H4 at largest optical 
depth is an artifact of the averaging procedure and should 
be ignored. 

vides the correct entropy jump. The actual run of the en- 
tropy in the optically thick layers is not very well matched: 
usually, a RHD model predicts a more rapid switching be- 
tween adiabatically and radiatively stratified layers. In M- 
type objects, a calibrated MLT model matches the actual 
thermal profile in the optically thick regions more closely. 
This property is likely related to the more gradual transi- 
tion between the two modes of energy transport. 

4.1. Spectroscopic effects 

In this section, we want to demonstrate which impact the 
differences between the thermal structur es of RHD and 
MLT models have on spectral properties. iMohantv et alJ 
l(2004l) used molecular bands of titanium-oxide and lines 
of neutral atomic alkalis to determine the effective tem- 
peratures and surface gravities of M-type PMS objects by 
comparing synthetic and observed spectra. The tempera- 
tures and gravities of the objects studied by Mohanty and 
collaborators fell into the regime considered here. The au- 
thors emphasized that the strengths of the investigated 
TiO band heads serve as excellent and important tem- 
perature indicator. Figures 1181 and 1191 show two prime 
spectral regions (wavelengths are given as wavelengths in 
air) considered in the analysis by Mohanty et al. . For our 
comparison, we picked case C3 where we found signifi- 
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Fig. 18. Comparison of synthetic spectra showing the 
triple-headed epsilon {E'^Ii - X^A) band of TiO (at 
8432, 8442, 8452 A), based on the hydrodynamical struc- 
ture C3 (black solid line) and a corresponding mixing- 
length model with omlt — 1-5 (grey solid line, red in 
color version) . Ti i lines are also see to absorb in the 8432 
epsilon subband at AA 8435.7, 8435.0 A, and bluewards of 
the triple band system at 8412.3 and 8426.5 A. Another 
strong Til line is also absorbing at 8468.4 A (not seen in 
this plot). 
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Fig. 19. Comparison of synthetic spectra showing a subor- 
dinate Nai doublet (at 8183.3, 8194.8 A) based on hydro- 
dynamical model C3 (black solid line) and a corresponding 
mixing-length model with omlt = 1-5 (grey solid line, red 
in color version). 



cant differences in the thermal structures between RHD 
and MLT models. Spectral synthesis calculations were per- 
formed with the PHOENIX code on the prescribed struc- 
tures at a spectral resolution of 0.01 A. In Figs. 1181 andll9l 
the spectral resolution has been degraded to ^ 30 000 sim- 
ilar to the one in the work of Mohanty et al. . 

We find that, accounting for the hydrodynamical struc- 
ture, yields systematically weaker TiO (and H2O not 
shown) bands by 0.18 and 0.025 dex respectively, while 
the pseudo-continuum appears unchanged. This is possi- 
ble because the strongest TiO bands are formed at two 
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dex lower optical depth than the opacity minima between 
those bands. The differences in the strength of the ep- 
silon subband heads in the synthetic spectra of Fig. [T51 
would correspond to a difference in Tgs of ^ 200 K when 
compared to an observed spectrum, in the sense that an 
analysis based on MLT models overestimates Teg. 

Atomic lines absorbing through TiO bands troughs 
such has the doublets of Ti at AA 8435.7, 8435.0 A, of Ki 
at AA 7664.9, 7699.0 A, and all other atomic lines formed 
bluewards of 0.7 /xm, would look wider and deeper in con- 
trast to the TiO pseudo-continuum, causing MLT mod- 
els to overestimate gravities. This is not the case of the 
Nal doublet at AA 8183.3, 8194.8 A shown in Fig. [Bl which 
is practically unaffected by this pseudo-continuum be- 
cause it forms between TiO and VO band heads from 
deeper photospheric layers. The same would be true of 
lines formed around the peak of the spectral distribution 
between 0.9 and 1.3 /im. Although of course these can be 
as well affected in analysis where TcS is determined from 
the TiO bands. 

5. Mixing by atmospheric overshoot 

For describing the mixing properties of the flow field in 
the overshooting layers of our models, we follow the ap- 
proximate procedure laid out by LAH. We describe the 
mixing in terms of a mass exchange frequency given by 

/pup 



(mcoi)(2;) 

^mass the upward directed component of the mass flux 



-^massC-^i ^1 ^) 



pv^ if Vz > 
otherwise. 



(3) 



where is the vertical component of the velocity (counted 
as positive if directed upwards), p the mass density, a;, y, 
z the spatial coordinates, and t the time, mcoi is the mass 
column density given by 



mcoi{x,y,z,t) 



dz' p(x, y, z', t). 



(4) 



(.) denotes the horizontal and temporal average over x, y, 
and t. The basic idea is to take the time scale over which 
the mass above a certain reference height is potentially 
exchanged by the flow as time scale over which material is 
mixed with fresh material stemming from the deeper ly- 
ing, convective layers. As we shall see the mass exchange 
frequency /cx exhibits an exponential height dependence. 
The mixing rate given by relation ^ is an approximation 
only. Depending on the way the mixing takes place in de- 
tail, the normalization of the mixing profile might change. 
However, the relative shape of the mixing rate - the expo- 
nential decline - is a robust feature, and in the following 
we shall characterize the mixing found in our models by 
the scale height of the exponential decline. 

We note that in Ludwis (200^ the mixing was de- 
scribed in terms of a mixing velocity 

\ mass/ 



(5) 



where p denotes the mass density. Wmix also shows an ex- 
ponential height dependence, and its rate of decline was 
given in the above paper. While at the level of our ap- 
proximation the description in terms of fmix is equivalent 
to the one in terms of /ex, the scale heights of the various 
declines can differ substantially. Differences become large 
in cases where the scale heights of Wmix or /cx are large, 
or the scale heights of p and mcoi differ noticeably. 

5.1. Subsonic filtering 

The atmospheric velocity field is a superposition of advec- 
tive motions and acoustic waves gen erated by convection 
in de eper layers (see also LAH, and iLudwig fc Nordliiiidl 
l200nl) . The wave motions contribute little if at all to the 
mixing due to their spatially coherent, oscillatory char- 
acter. The overshooting, convective motions tend to de- 
cay with distance to the Schwarzschild stability boundary, 
while the wave amplitudes tend to increase with height 
due to the sharp decrease in mass density. This leads to 
the situation that beyond a certain height the atmospheric 
velocity field becomes dominated by wave motions. In or- 
der to get a reliable estimate of the mixing, it is therefore 
necessary to remove the wave contributions to the velocity 
field before evaluating the mass flux F^^^. 

We removed the wave contributions by subsonic filter- 
ing - a technique developed in the context of solar obser- 
vations for cleaning images f rom "noise" stem ming from 
the solar 5 minute oscillations l|Title et al.ll989i) . FigureEHl 
schematically illustrates this filtering technique. In short, 
one considers a time sequence of images and removes fea- 
tures with horizontal phase speeds Wphase greater than a 
prescribed threshold. This is achieved by Fourier filter- 
ing of spatial-temporal data in the k-tu domain. For every 
depth layer in our data cubes we performed a 3D Fourier 
analysis (one temporal, two spatial dimensions) of the ver- 
tical mass flux retaining only contributions below a preset 
phase speed threshold. In practice, acoustic and convec- 
tive contributions are not as cleanly separated as shown in 
the Fig. 1201 and one must find the right balance between 
removing as much acoustic components as possible while 
retaining as much as possible convective contributions. We 
always studied a sequence of phase speed thresholds in or- 
der to judge the success of the procedure. 

5.2. IVIixing in the solar atmosphere 

For reference we begin with a discussion of the mixing 
in the solar atmosphere. Figure 1211 shows /cx in our so- 
lar model for various degrees of subsonic filtering. It is 
clearly visible that the subsonic filtering has the strongest 
impact on /cx in the uppermost atmospheric layers. As 
hinted above, /cx exhibits an exponential decline with 
height (logP cc z) afte r appropriate s ubson ic filtering 
which was put forward bv lFrevtag et al.l l)l996j) as generic 
feature of convective overshoot. Figure TIT] further shows 
that too low a velocity threshold removes also convective 
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Fig. 20. Schematic illustration of subsonic filtering in the 
k-iij domain: only components in wavenumber-frequency 
domain below a prescribed phase speed Wphasc are retained. 
They preferentially belong to convective motions. 



features, as visible by the reduction of the velocity in the 
deeper, convection dominated layers for the case of a phase 
speed threshold of 1.5kms~^. Quantitatively, for the Sun 
we find a scale height of /ex in terms of the local pres- 
sure scale height of iJfcx = 2.4 Hp with an uncertainty'^ of 
about 10%. 

Apart from purely numerical findings, an exponen- 
tial decline of /ex is also motivated from semi-analytical 
considerations of the behavior of linear convective modes 
i Frevtag et 31111996(1 . In Fig. EH we plotted the /ex-profile 
of a linear convective eigenmode with horizontal wave- 
length of Amodo = 5.0 Mm. We used the temporally and 
horizontally averaged hydrodynamical structure as back- 
ground on which we solved the linearized hydrodynamical 
equations. The absolute amplitude of the mode has been 
scaled to match /cx of the hydrodynamical model leaving 
the shape of the mode's /ex-profile intact. The exchange 
frequency of the mode exhibits an exponential "leakage" 
into the formally convective stable layers. Generally, the 
rate of decline depends on Amodo, being faster for modes 
of shorter horizontal wavelengths. The mode with 5.0 Mm 
wavelength was chosen since it provided a good overall fit 
to the decline of /ex in the hydrodynamical model. The 
wavelength of this mode is significantly larger than the 
horizontal scale of the dominant convective structures on 
the Sun - the granules with typical sizes of around 1.2 Mm. 
This might be related to the assumption of adiabaticity in 
the mode calculations, which is not a good approximation 
in the solar photosphere, or to the fact that a convective 
mode is a non-stationary solution of the hydrodynamical 
equations. 



^ The uncertainty is not meant in a statistical sense but re- 
flects the precision with which we can read off the slope from 
the plots of /ex. The alert reader mig ht suspect a conn ection 
to the "chi-by-eye" technique, consult IPress et alJ ilQQST) for a 
discussion of the immediate consequences. 
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Fig. 21. Mass exchange frequency /ex in solar model S as 
a function of gas pressure: the unfiltered data (thick solid 
line) were subsonically filtered retaining only features with 
phase speeds Uphasc < 12, 6, 3, and 1.5kms~^ (thin solid 
lines from top to bottom) . The approximate location of the 
Schwarzschild boundary of convective stability is indicated 
by the dotted vertical line, /ex of a convective eigenmodes 
with horizontal wavelength Amodc = 5.0 Mm = 33i/p"''^ 
by the dashed line. The dash-dotted line is a fit depicting 
the decline of /ex- It is labeled by the scale height of the 
decline in units of the local pressure scale height. The 
pressure is given in units of the pressure at Rosseland 
optical depth unity Psurf • The plateau at lowest pressures 
is an artifact of the upper boundary condition applied in 
the RHD model. 



From the rather large wavelength of the best fitting 
mode one might argue that 5.0 Mm is close to the geomet- 
rical size of the computational box of model S (6.0 Mm), 
and actually the box size sets the rate of decline of /ex- 
We verified that a solar model of about twice the hori- 
zontal size gives the same rate of decline as model S. The 
box size of model S is sufficient to allow the build-up of 
all convective structures contributing significantly to the 
overshooting velocity field in the deep photosphere. The 
box sizes of the M-type models are allowing the presence 
of a similar number of convective cells as the solar model. 
Thus, we expect that also our M-type models capture the 
relevant convective structures controlling the overshooting 
motions. 

5.3. Mixing in M-type atmospheres 

In Fig. 122 and we show the vertical distribution of 
the mass-exchange rates /ex for two of our M-type mod- 
els. With decreasing logg and increasing Teff, the zone 
of convective instability extends further and further into 
the optically thin atmosphere, leaving little room for over- 
shoot in the models C3 and H4. Reading off an exponential 
decline rate is very uncertain in these models. However, 
from them we find a slow decline with -fffex ~ 18 Hp in C3 
and w 28 Hp in H4. Models C3 and C4 leave more room 
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for overshooting, allowing a more precise determination 
of the exponential decline. We estimate the involved un- 
certainty to about 20%. We find fffex = 0.5 -ffp in model 
C5 and 3.2 Hp in C4. Note, that in the model C5 with 
highest gravity and steepest decline of /ex the exponential 
behavior does not set in immediately at the boundary of 
convective stability 

Qualitatively, in terms of the rate of decline, overshoot- 
ing is less pronounced in models of higher gravity. This is 
in part due to the fact that buoyancy forces scale pro- 
portional to gravity, making buoyancy more effective in 
confining the convective motions to the formally unstable 
regions. At lower gravity, mixing - despite the increas- 
ing geometrical scales - is more rapid, not only due to 
the slower decline of the mixing rate but also due to the 
higher convective velocities. 

Silicate cloud formation is one of the most important 
aspect of the modeling of late- type M and brown dwarfs. 
The formation of clouds is understood as a compromise 
between condensation, sedimentation and advection (tur- 
bulent overshooting mixing) time scales which determine 
the extension, location in the thermal atmospheric struc- 
ture, and composition of the cloud deck. To represent the 
correct distribution with height of the mixing time scale, 
investi gators hav e experimented with various descriptions 
for f.,^: I AUard et al. (2003) used a parabolic function with 
opening set by the innermost and outermost convective 
layers, and normalized at the convect ive velocity maxi- 
mum, while IXckerman fc MarlevI l)200l|) and lCoooer et all 
l|200.'^ preferred a constant distribution throughout the 
atmosphere, set to the value associated with the maxi- 
mum of the convective velocity. A modeling with convec- 
tive modes should give a more physical description which 
could be implemented in ID model atmospheres. Trying 
to match the mixing profiles in the overshooting regions 
with convective modes, however, worked only partially so 
far. For the lower gravity models the fits were not satis- 
factory. This might be related to the situation that con- 
vection reaches high up, and we do not actually see the 
asymptotic exponential tail of the mixing profile. We ori- 
ented the horizontal wavelength of the linear modes at the 
largest sizes of structures the computational box could ac- 
commodate in the respective models. Despite the present 
shortcoming we are optimistic that one can add refine- 
ments to the mode-modelization that would allow to sat- 
isfactorily match the RHD results. 

If the mixing trends observed in our models hold for 
cooler objects, these go in the direction of making clouds 
thicker or more extended into higher atmospheric layers 
with decreasing gravities and increasing T^g . However, de- 
creasing pressure will work in the opposite direction, mak- 
ing it harder for grains to form. Detailed calculations will 
be presented in a subsequent publication. Nevertheless, 
we expect that clouds will be more extended for young 
objects than for older ones of same T^g, and that these 
will remain dusty at lower Toff and la ter spectral types, 
i.e, b elow spectral class T4 or 1400 K (|Golimowski et all 
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5.4. Cloud cover disruption in early T-type brown 
dwarfs 

iBurgasser et al.l ()2002f) have found a resurgence of molec- 
ular spectral features such as FeH bands in the spectra 
of early T-type brown dwarfs. This is interpreted as a 
spectral signature of the onset of cloud cover disruption. 
Indeed, these spectral features of refractory species can 
only be seen if the atmosphere is transparent enough to 
observe flux emerging from below the cloud forming layers. 
This is possible if holes in the cloud deck are occurring. 

Dust does not form in the models studied in this work. 
However, here we want to speculate how the cloud pattern 
might look when one expects a disrupted cloud layer like 
in early T-dwarfs. The cloud deck is shaped by convec- 
tive overshooting which mixes up refractory material into 
the grain condensing part of the atmosphere - below gas 
temperatures of 2000 K, in early T-dwarfs perhaps over 
one pressure scale height above the convectively unstable 
layers. This far above the convection zone, the horizontal 
and vertical motions are not correlated in the same way as 
in the strongly convective layer where the flow forms cell- 
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like patterns. Structures larger than the granular scale in 
combination with waves dominate the velocity field. The 
typical granular flow pattern is "washed out" from the 
flow higher up. Hence, we do not expect that the cloud 
deck is fragmented on a spatial scale given by the granu- 
lar scale, but likely on a larger scale. This consideration 
refers to effects of convection. It is of course well possi- 
ble that the actual cloud pattern is rather shaped by the 
global wind circulation expected to be present in rotating 
brown dwarfs or planets. 

6. Final remarks 

We have seen that mixing- length theory provides a reason- 
ably realistic picture of the convective energy transport in 
M-type atmospheres, even considering that a substantial 
part of the optically thin atmosphere is affected by con- 
vection. However quantitatively, temperature errors of up 
to « 250 K (or 9 %) are possible if one (unluckily) picked a 
value of unity for the mixing- length parameter (aMLx) en- 
tering MLT. The efficiency of the convective energy trans- 
port measured in terms of an effective mixing-length pa- 
rameter is rather high in M-type atmospheres. Choosing 
a larger c^mlt helps but is not sufficient to describe the 
thermal structure of M-type atmospheres if one wishes 
to attain a high level of accuracy. MLT does not provide 
the precise scaling of the convective transport efficiency 
with stellar parameters and optical depth. To get a better 
quantitative description one might try to calibrate besides 
Q!MLT the "internal" parameters of MLT with radiation- 
hydrodynamics models. M-type atmospheres appear par- 
ticularly well suited for this undertaking since convection 
takes place under optically thick and thin conditions. 

Our results have shown that convective overshooting 
mixes the layers of M-type atmospheres which are formally 
(according to the Schwarzschild criterion) stable against 
convection much more strongly towards higher effective 
temperatures and lower gravities. This should make young 
late-type M-dwarfs and brown dwarfs even more cloudy 
than older disk objects of the same Tcs- We think that a 
detailed modeling of such cooler atmospheres, especially 
around 1400 K should deliver important clues about the 
interesting question of the cloud cover, and should help 
to understand the break in colors and spectral type vs. 
effective tem perature relation obse rved for these objects 
( KnapD et ah 2004: Golimowski et al. 200 4) . 

To aid the modeling of the spectral properties of cloudy 
objects, we further think it should be worthwhile to im- 
prove the model of convective modes which where mostly 
used for demonstration purposes in this work. Our current 
mode model did not perform sufficiently to be up to the 
task but a number of refinements can be brought to this 
model and tested. It could be calibrated with RHD models 
in a similar fashion like MLT. At a higher ambition level 
one might even contemplate to combine ideas into a single 
description for the energy transport and effects related to 
overshooting. We are aware that many attempts have been 
made to improve or even completely replace MLT since it 



has been introduced into astrophysics in the early 1950s 
by Bohm-Vitense - with mixed success. Our goal would 
not be to formulate a new convection theory but rather 
a parameterized model like MLT which is flexible enough 
to fully fit RHD results, and contains sufficient physics to 
allow a robust inter- and extrapolation in a wide range 
of stellar parameters. The availability of detailed RHD 
models appear essential to identify the necessary building 
blocks. 

The main uncertainty affecting our present results is 
related to the approximate treatment of the wavelength- 
dependence of opacity in the radiative transfer which 
was optimized for an atmosphere at Tefi=2900K and 
logSf=5.3. So, one of the first issues to be addressed in 
future work is the improvement of this approach. Work 
is under way for a refined implementation in a new 
3D radiation-hydrodynamics code (named CO^BOLD, 
iFrevtag et all l2002t IWedemever et al.l l2004j) . We finally 
emphasize that our results apply to atmospheres of solar 
metalicity. We expect marked differences f or metal-poor 
atmospheres (see, e.af.. lAsnlund et aDll999l) . 
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Appendix A: Computation of the Peclet number 

The Peclet number Pe measures the relative importance 
between conductive (here by radiation) and advective heat 
transport 

^rad 



Pe 



adv 



(A.l) 



trad is a radiative relaxation time, and tadv a characteristic 
time over which the temperature of moving gas elements 
changes due to adiabatic compression or expansion. In the 
present context, we employ a mixing-length picture and 
evaluate the radiative relaxation time t^ad with the MLT 
formula 



^rad 



(A.2) 



Cp denotes the specific heat at constant pressure, A = 
Q^MLT^^p the mixing-length, a Stefan-Boltzmann's con- 
stant, X opacity, T temperature, p mass density, and 
the optical thickness of a convective element defined as 

Tc = XPA. (A.3) 
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fs — 16 and f4 — 2 are dimensionless constants set to 
values assumed in the MLT formulation of Mihalas (see 
iLudwig et alJll999l) . We further assume a mixing-length 
parameter omlt = 2.5 which is a reasonable value for the 
M-type atmospheres under consideration (see Fig. [T^ . 

Similarly, we estimate tadv as the time interval over 
which a vertically moving gas element has build up a sub- 
stantial temperature difference according to 

iadv = — A.4) 

where V is the logarithmic temperature derivative of the 
thermal profile with respect to pressure, Vad the corre- 
sponding adiabatic value, and Vc a convective velocity we 
set to a typical atmospheric value (300 m s^^ for model 
C5 and 600m s"^ for model C3, see Fig.[l7|). Note, that 
in this paper we argue taking recourse to ratios of Pe 
only. This makes the precise choice of arbitrary or little 
constrained parameters less critical. 

References 

Ackerman, A. S. & Marley, M. S. 2001, ApJ, 556, 872 
Allard, P., Guillot, T., Ludwig, H.-G., et al. 2003, in lAU 

Symposium No. 211: Brown Dwarfs, 325 
Asplund, M., Nordlund, A., Trampedach, R., & Stein, 

R. F. 1999, A&A, 346, L17 
Baraffe, L, Chabrier, G., Allard, P., & Hauschildt, P. 2002, 

A&A, 382, 563 
Bohm-Vitense, E. 1958, Zs. f. Ap., 46, 108 
Burgasser, A. J., Marley, M. S., Ackerman, A. S., et al. 

2002, ApJL, 571, L151 
Cattaneo, P., Brumell, N., Toomre, J., Malagoli, A., & 

Hurlburt, N. 1991, Astrophys. J., 370 
Chan, K. & Sofia, S. 1987, Science, 235, 465 
Chan, K. & Sofia, S. 1989, Astrophys. J., 336, 1022 
Cooper, C. S., Sudarsky, D., Milsom, J. A., Lunine, J. L, 

& Burrows, A. 2003, ApJ, 586, 1320 
Perguson, J. W., Alexander, D. R., Allard, P., et al. 2005, 

Astrophys. J., 623, 585 
Freytag, B., Ludwig, H.-G., & Steffen, M. 1996, A&A, 313, 

397 

Freytag, B., Steffen, M., & Dorch, B. 2002, AN, 323, 213 
Gadun, A. & Pikalov, K. 1996, Solar Phys., 166, 569 
Golimowski, D., Leggett, S., Marley, M., et al. 2004, AJ, 
127, 3516 

Hauschildt, P., Allard, P., & Baron, E. 1999, Astrophys. 
J., 512, 377 

Helling, C, Klein, R., Woitke, P., Nowak, U., & Sedlmayr, 

E. 2004, AA, 423, 657 
Kim, Y.-C, Pox, P., Sofia, S., & Demarque, P. 1995, 

Astrophys. J., 442, 422 
Knapp, G., Leggett, S., Pan, X., et al. 2004, AJ, 127, 3553 
Ludwig, H.-G. 1992, Ph.D. Thesis (University of Kiel) 
Ludwig, H.-G. 2003, in Proceedings of the 210*'' lAU 

Symposium: Modelling of Stellar Atmosperes, ed. 

N. Piskunov, W. Weiss, & D. Gray (Astronomical 

Society of the Pacific), 113 



Ludwig, H.-G., Allard, P., & Hauschildt, P. 2002, A&A, 
395, 99 

Ludwig, H.-G., Freytag, B., & Steffen, M. 1999, A&A, 346, 
111 

Ludwig, H.-G., Jordan, S., & Steffen, M. 1994, A&A, 284, 
105 

Ludwig, H.-G. & Nordlund, A. 2000, in Pacific Rim 
Conference on Stellar Astrophysics, ed. K. Cheng, 
H. Chau, & K. Chan (Kluver Academic Publishers), 
37-44 

Mihalas, D. 1978, Stellar Atmospheres (Freeman and 
Company) 

Mohanty, S., Basri, G., Jayawardhana, R., et al. 2004, 

Astrophys. J., 609, 854 
Nordlund, A. 1982, A&A, 107, 1 
Nordlund, A. & Dravins, D. 1990, AA, 228, 155 
Press, W., Teukolsky, S., Vetterling, W., & Plannery, 

B. 1992, Numerical Recipes in Fortran (Cambridge 

University Press) 
Robinson, F., Demarque, P., Li, L., et al. 2003, MNRAS, 

340, 923 

Robinson, P., Demarque, P., Li, L., et al. 2004, MNRAS, 
347, 1208 

Singh, H. & Chan, K. 1993, A&A, 279, 107 

Steffen, M. 1993, in Inside the stars, ed. W. Weiss & 

A. Baghn (ASP Conference Series, vol. 40), 300 
Steffen, M., Ludwig, H.-G., & Freytag, B. 1995, A&A, 300, 

473 

Steffen, M., Ludwig, H.-G., & KriiB, A. 1989, A&A, 213, 
371 

Stein, R. P. & Nordlund, A. 1998, Astrophys. J., 499, 914 
Steiner, O., Grossmann-Doerth, U., Knoelker, M., & 

Schiissler, M. 1998, ApJ, 495, 468 
Title, A. M., Tarbell, T. D., Topka, K. P., Ferguson, S. H., 

& Shine, R. A. 1989, Astrophys. J., 336, 475 
Tsuji, T. 2005, ApJ, 621, 1033 

Vogler, A., Bruls, J., & Schiissler, M. 2004, AA, 421, 741 
Vogler, A. & Schiissler, M. 2003, AN, 324, 399 
Wedemeyer, S., Freytag, B., Steffen, M., Ludwig, H.-G., 
& Holweger, H. 2004, A&A, 414, 1121 



